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1 Introduction and discussions 

Recent years have witnessed major breakthroughs in computing and understanding scat¬ 
tering amplitudes of gauge theory and gravity (see, e.g., [1] for reviews). At many stages 
of the development, the AT = 4 supersymmetric Yang-Mills theory (SYM) proved to be an 
extremely fruitful testing ground. Many novel ideas, such as twistor string theory [2], dual 
superconformal symmetry [3], Grassmannian formulation [4], on-shell diagram representa¬ 
tion [5] and amplituhedron [6, 7], are realized in their simplest forms in AT = 4 SYM and 
then generalized to less symmetric theories. 

The three dimensional AT = 6 supersymmetric Chern-Simons-matter theory [8-14], of¬ 
ten referred to as the ABJM theory, is a close sibling of the AT = 4 SYM in many respects. 
For instance, the scattering amplitudes of the ABJM theory exhibit dual superconformal 
symmetry [15-18] and admits Grassmannian [19], twistor string [20-22], and on-shell dia¬ 
gram formulations [23-25]. Despite these parallel successes, the study of ABJM amplitudes 
fall short of those of AT = 4 SYM in many respects. One of the most pressing problem is 
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the lack of a “momentum twistor” formulation [26-28] in which the dual super conformal 
symmetry would become manifest (see [29] for a recent attempt in this direction). Closely 
related to momentum twistors are the dual superconformal ill-invariants [3], which serves 
as building blocks for an explicit formula for all tree amplitudes [30] and a starting point 
for the construction of amplituhedron. 

In principle, all ingredients to compute the ABJM tree amplitudes are available in the 
literature. The Grassmannian integral [19], supplemented by the contour prescription from 
on-shell diagrams [23-25], will produce the amplitudes. A mundane, yet seemingly unavoid¬ 
able, problem is that each BCFW bridge in the on-shell diagram introduces a quadratic 
equation in the integration variables. The solutions to quadratic equations generically con¬ 
tain square-roots, which must cancel out when summed over all solutions and produce a 
rational function of kinematic variables. Mainly for this technical reason, explicit results 
for ABJM tree amplitudes to date are limited to 4- and 6 -point amplitudes [15, 18, 19] 
which are free from square-roots due to limited kinematics, and a partial result for 8 -point 
amplitude [18] without manifest supersymmetry. 

The goal of this paper is to take a few steps toward the computation of all ABJM tree 
amplitudes. Our two main results are a complete evaluation of the supersymmetric 8 -point 
amplitude and a derivation of a double soft theorem valid for all tree amplitudes. 

In evaluating the Grassmannian integral for 8 -point or higher amplitudes, we find it 
convenient to use a new gauge, which we call “u-gauge”. The u-gauge is inspired by the 
isomorphism between the orthogonal Grassmannian and the pure spinor geometries; both 
of them admit the SO{2k)/U(k) coset description. A particular set of coordinates of the 
coset space introduced in [31] trivially solves the orthogonality constraint and can be easily 
generalized to arbitrary k. This fact makes the u-gauge, at least in some contexts, more 
convenient than conventional gauges involving Euler angle coordinates. 

Although the u-gauge do not circumvent the square-root problem mentioned above, 
the quadratic equations in the u-gauge tend to be simpler, which allow us to combine all 
residues in the contour integral. For the 8 -point amplitude, the integral is effectively one¬ 
dimensional. We can express the denominators of the amplitude in terms of the standard 
cross-ratios among solutions to quadratic equations. It is easy to see that the cross-ratios 
can in turn written in terms of the coefficients of the quadratic equations, thereby avoiding 
the need to solve the equations explicitly. 

The final result for the 8 -point amplitude takes the form, 



( 1 . 1 ) 


The (1 + 7 r) factor accounts for the sum over two disjoint branches of the orthogonal 


Grassmannian. The two rational functions in the big parenthesis corresponds to the two 


on-shell diagrams contributing to the 8 -point amplitude. The numerators F’(l), .F(3) as 
well as the A t j factors in the denominators are polynomials in kinematic variables. The 
on-shell diagrams suggest that A 12 , A 14 , A 32 , A 34 should be proportional to physical poles 
of the amplitude whereas A 13 = A 31 should be spurious. We confirm the expectation 
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by explicitly proving that A^ for the physical poles are proportional to pf.i rn factors for 
adjacent particles. 

In the second half of this paper, we consider the double-soft limit of ABJM ampli¬ 
tudes. Soft limits of scattering amplitudes in gauge and gravity are well known to exhibit 
universal behavior and have bearing on gauge symmetries and spontaneously broken global 
symmetries. Our motivation to study the soft limit is more modest. As we make progress 
in computing higher point amplitudes, we wish to use the soft theorem to test the consis¬ 
tency of the methods we use. Our derivation of the soft theorem will closely follow that 
of ref. [32], where a similar double-soft theorem was derived for three dimensional super¬ 
gravity theories. We show that the (2k + 2)-point amplitude A 2 k +2 reduces to the 2A;-point 
amplitudes A 2 k with universal leading and sub-leading soft factors in the double soft limit, 

-^ 2 fc+ 2 (l 1 2 , • • • ,2k, e 2 p2k+i, e 2 P 2 fc+ 2)| e _> 0 = + £ S (1) ) -4.2fc(l> 2 , • • • , 2k) . ( 1 . 2 ) 

As in [32], the proof of the soft theorem is based on the BCFW [33, 34] recursion relation of 
the ABJM theory [18]. We confirm that the universal soft factors respect all the symmetries 
of the ABJM amplitudes for all k. For 6 -point amplitude to the sub-leading order, and 
for 8 -point amplitude to the leading order, we explicitly take the soft limit of the known 
amplitude and verify that the soft theorem holds. 

Although we still have explicit form of tree amplitudes only up to 8 -point, we expect 
that the findings in the present paper, such as the u-gauge, cross-ratios among different 
poles in the contour integral, and the double soft theorem, will lay the groundwork for 
a complete construction of all tree amplitudes of ABJM theory in terms of momentum 
twistors and/or dual superconformal ^-invariants. 

This paper is organized as follows. In section 2, we give a short review of the general 
structure of the ABJM tree amplitudes and the Grassmannian integral. Then we introduce 
the u-gauge and compare it with other well-known gauges. In section 3, we use the u-gauge 
to compute some tree amplitudes. After reproducing the 4-point and 6 -point amplitudes, 
we present the details of how to evaluate the 8 -point amplitude. In section 4, we propose 
the double soft theorem of the ABJM amplitudes and prove it using the BCFW recursion 
relation. We take the double soft limit of the 6 - and 8 -point amplitudes, and verify explicitly 
that the theorem holds. 

2 Grasssmannian integral in the U-gauge 
2.1 Preliminaries 

Here we briefly review the 3-dimensional spinor helicity formalism [15] and introduce the 
Grassmannian integral formula for tree level amplitudes of planar ABJM theory [19]. 

Each spinor in three dimensions transforms under SL(2, M), and a null momentum can 
be written in the bi-spinor form 

p afl = Pfl (c7^ = \ a \P. (2.1) 
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Our convention for spinors and gamma matrices are such that p a & is real for real p^, and 
A" is real (purely imaginary) for outgoing (incoming) particles. The spinors are contracted 
as (ij) = XfX j a . We normalize the norm of vectors such that p'X = (pi + Pj) 2 = {ij) 2 when 
both A i and A j are real. 

The on-shell superheld notation for ABJM amplitudes is built on three fermionic coor¬ 
dinates r] 1 , in addition to X a , which transform as 3 under the U(3) subgroup of the SO( 6 ) 
R-symmetry group. The particle/anti-particle superfields take the form 

$> = <t> 4 + r/'ipr + \e IJK r] I p J (l) K + 4, (2.2) 

2 b 

$ = + rffa + ]-e IJK ri I ri J/ ip K + \ti j W r] J r] K (2.3) 

2 o 

A collective notation A = (A; 77 ) will be used when appropriate. The fact that (2.1) is 
invariant under A —> — A, while the wave-functions of fermions pick up a minus sign, implies 
the so-called “A-parity” of the super-amplitudes. 

Al(Ai,..., —A i ,..., A 2k) = (—l)M(Ai,..., A i ,..., A 2 fc). (2.4) 


The super-conformal generators of the superconformal symmetry come in three types: 

d . . d 2 


A 


dA’ 


AA, 


9AM’ 


(2.5) 


For the second type, we will use the notation 

p aP = x a x ^ q od = ( 2 . 6 ) 

The super-momentum conservation is denoted as delta functions by 

2k 2k 

5 3 (P)5 6 (Q) with P = J^pf, ( 2 - 7 ) 

i =1 z=l 

The Grassmannian integral formula for the tree level amplitudes of planar ABJM 
theory, first proposed in [19], is 


/• d kx2k C 5 k ( k +i)/Z(C-C T )5 2k \ 3k {C ■ A) 

C2k{ ’ ~ J vol[GL(k)] M x M 2 ---M k _ x M k ' (2-8) 

The integration variable C is a (k x 2k) matrix. The dot products denote ( C ■ C T ) mn = 
Cm.iC n i , (C • A) m = C m iAi. The consecutive minor Mj is defined by 

Mi = e mi • • • C mk (i +k _iy (2-9) 

It was shown in [19] that this formula satisfies the same cyclic symmetry and superconformal 
symmetry as the tree-level (2k)-point amplitude. Yangian invariance of the formula was 
first argued in [19] and explicitly proved later in [35]. 
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The integral (2.8) should be considered as a contour integral on the moduli space of 
rank k, ( k x 2k) matrices C with the constraint C -C T = 0 and the equivalence relation C ~ 
gC (g E GL(A;)). This moduli space is known as the orthogonal Grassmannian OG(fc, 2k). 
The dimension of OG (k, 2k) is determined by the aforementioned two conditions: 

dim c [OG(fc, 2k)] = 2 k 2 - k 2 - k i k ±Jl = ill. (2.10) 

Integrating out the bosonic delta function d 2k (C ■ A) leaves the momentum conserving 
delta function and a contour integral over {k — 2 )(k — 3)/2 variables. The geometry and 
combinatorics behind the Grassmannian integral for all tree amplitudes, as well as some 
loop amplitudes, have been elucidated in [23-25]. On the other hand, explicit computation 
of amplitudes has never proceeded beyond 8-point [18]. 

2.2 U-gauge 

We take a real slice of the complex orthogonal Grassmannian with the split signature, 
where the “metric” in the particle basis is 

g = diag(—, +,—,+,•••)• (2.11) 

In this basis, the momenta and their spinor variables are related by 

pf = (-l)'A?Af, Af€R, (P,+Pi) 2 = (-1) ,+I (ij) 2 - (2.12) 

We find it convenient to switch between the particle basis and the light-cone basis: 

2k k 

ds — ^ ( ( 1 ) (dX{) — ^ ( dw m dv — ^ 2 m T 3'2m— I; V — m % 2 m— l) • (2.13) 

i= 1 m=l 

In the light-cone basis, a GL(k)n subgroup of the 0(k,k) symmetry group remains man¬ 
ifest. We will use a notation with covariance under SL(k) C GL(k)n and adopt the 
summation convention. The invariant tensor of SL(k,M) will be denoted by 

The spinor-helicity variable for particles are denoted by Xf (i = 1,, 2k, a = 1, 2). 
We use the same letters w, v for the light-cone combinations of the spinor variables: 

w m = ^2m + ^2m—l > ' L ’ ma = ^2m ~ ^2m-l ■ (2-14) 

The scalar product of two spinors are defined in a usual manner. 

(AiA 2 ) = e Q/3 A ?\% . (2.15) 

To avoid confusion, we reserve the shorthand notation (12 ) = (AiA 2 ) exclusively for the 
particle basis. In the light-cone basis, we will use (w m w n ), (w m v n ) and so on. The overall 
momentum conservation is written as 

2k 

P 0/3 = E(-l)'Af A? = j (<^ + v&v ma ) = 0 • ( 2 -!6) 

i= 1 
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The light-cone components of the fermionic coordinates gj are denoted by 

C, = vim + vL-i > = ? ?L - vL-i ■ 

The supermomentum components are rewritten as 


■2k 


1 


Q aI = £(-l)W = 2 ' 


2=1 


(2.17) 


(2.18) 


In summary, the metric and the kinetic variables in the light-cone basis take the form 


9 = 


1/0 5 m 7 
2 S m n 0 


, A = 


W r , 


V=\ 0n 


(2.19) 


The light-cone form of the C-matrix before a gauge fixing is 

C=(t„ n 


'a a an 


( 2 . 20 ) 


A priori, the GL(k)i index a is not correlated with the light-cone index n. We choose to 
fix the gauge by locking GL(k)L and GL(k)n: 


C — | San'j t \ Umn ^ j ^ mn — (2 )m 

The orthogonality condition implies that u mn is anti-symmetric: 

5(C ■ g ■ C T ) = 5 ( \{Umn + u nm j) ■ 


( 2 . 21 ) 


( 2 . 22 ) 


Since the decomposition of u m n into the symmetric and anti-symmetric parts is a linear 
operation, the delta-function does not produce any u-dependent Jacobian factor. 

We will call this gauge fixing the “u-gauge”. This gauge was inspired by the fact that 
the orthogonal Grassmannian and the pure spinor admit the same SO(2k)/U(k) coset 
description and that the u mn coordinates were used in ref. [31] to solve the non-linear 
constraints of the pure spinors in order to construct higher dimensional twistor transforms. 

It is well known that C ■ A = 0 and C ■ g ■ C T = 0 implies the overall momentum 
conservation. In the light-cone gauge, C ■ A = 0 is written as 


tUm T 'U’mn'V — 0 • 

This equation admits a particular form of SL(k, R)-invariant solution for all k: 

* 2 (WmW n ) 1 (W p V P ) 

u = ~^{w m Wn) , R = --- ■ 


(■ WpVP) 


R 


(2.23) 


(2.24) 


To verify that (2.24) is indeed a solution to (2.23), it suffices to use the Schouten identity 


(w m w n )v pa + (w n v p )w “ + ( v p w m )w° = 0, 


(2.25) 


and the momentum conservation (2.16). 
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The light-cone basis before the gauge fixing respects the symmetry exchanging w n and 
v n . Thus it is natural to consider the “dual u-gauge” in which the roles of w n and v n are 
reversed: 

C = (x a n I u an ) ->• ( u mn | 5 m n ) . (2.26) 

In the dual u-gauge, the C ■ A = 0 condition reads 

v m + u mn w n = 0 , (2.27) 

which admits a particular solution, 

uT = ~\{v m v n ) . (2.28) 

The dual u-gauge will be useful in a later discussion on the A-parity for odd k. Using the 
energy momentum conservation and Schouten identity, one can show that 

< p u7 = -i( w m v n ), u* mp uT P = -2 . (2.29) 

For k = 2 and k = 3, (2.24) is the unique solution to (2.23). For higher k, there is a 
(k — 2){k — 3)/2-dimensional solution space containing (2.24). For instance, for k = 4, the 
general solution can be parametrized by 

'U-mniz) — Umn T mnrs • 

The general solution for k = 5 is 

z p 

'Umn(z') — ^mn T -^-tmnprs'U'* ■ 

The “vector” z p appears to have five components, but only three 
due to the equivalence relation, 


(2.30) 


(2.31) 

of them are independent 


z p + c a v- 


,pa 


(2.32) 


which follows from the fact that rt” s oc (v m v n ) and the Schouten identity. Along the same 
line of reasoning, we can write the general solution for k > 4 as 


u r , 


U mn + 


1 


2(k — 4)! 


Pl-Pfc —4, 


c 7rmpi ••’Pk—4'rs'U'* 


rs 


(2.33) 


In the (k — 2)-dimensional space surviving the quotient z p ~ z p + c a v pa , the tensor z pi '" Pk ~ 4 
spans a (k — 4)-plane. The effective number of components for z Pl "' Pk ~ 4 is, as expected, 


(k — 2\ _ fk - 2\ _ (k -2){k- 3) 

U-V “ V 2 ) - 2 


(2.34) 
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Broader class of u-gauges 

Most features of the u-gauge survives permutations of columns. For instance, instead of 
the alternating signature metric (2.11), we may take 

9 = dia g( r>''' > +>''' > + ) > ( 2 - 35 ) 

k k 

and define light-cone coordinates by 

w m = (-l) fc_1 (x m+ fc + x m ) , v m = (-1 ) fc_1 (x m+fc - x m ) . (2.36) 

One of the useful features of all u-type gauges is that the Jacobian factors arising from the 
computation of amplitudes are always powers of the R factor defined in (2.24), although 
the numerical value of R does depend on the particular gauge. 

There are (2 k)\/(k\) 2 different ways to distribute (—1) and (+1) in the diagonal entries 
of the metric. An overall flip of the signs is irrelevant, so there are (2k — 1 )\/(k\(k — 1)!) 
inequivalent metrics. Given a fixed metric, there are k\ inequivalent ways to pair the 
coordinates to define light-cone coordinates. To sum up, the number of different u-type 
gauges is (2k — 1 )\/(k — 1)!. 

Among all possibilities, we will mostly focus on the two choices we mentioned explicitly 
above. Both of them generalizes to arbitrary k straightforwardly. The alternating signature 
gauge defined (2.11), (2.13) is the only choice which respects the cyclic symmetry. For this 
reason we will call this gauge “u-cyclic gauge”. As we will see later, the other gauge defined 
by (2.35), (2.36) is convenient when we examine the factorization of A 2 k into two copies of 
Ak+i when k is odd. We will call this choice “u-factorization gauge”. 

Lambda-parity in the u-gauge 

Let us examine how the lambda parity is reflected in the u-gauge. We will show that, for 
odd k, the lambda parity induces the exchange, 

w m o v m . (2.37) 

For notational convenience, we will work in the u-factorization gauge, but the same argu¬ 
ments hold in all u-type gauges. 

With a usual gauge fixing in the particle basis [18, 19], the C-matrix is given by 

C=(l\6), OO t = I. (2.38) 

In the light-cone basis, the C-matrix translates to 

<7= (l-O | J + o) . (2.39) 

Note that the following identities hold for odd-dimensional orthogonal matrices: 

det(0) = ±l det(/=FO) = 0. (2.40) 
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When det(O) = — 1, a GL(k) gauge transformation gives rise to 

(l\ (/-0)- 1 (/ + 0)) = (j| Uj . (2.41) 

This establishes the relation between the u-gauge and the usual gauge in the particle basis. 

The other branch with det(O) = +1 is related to the det(O) = —1 branch by 
det(— O) = —det(O). The sign flip has the same effect as flipping the signs of all 
for m = k + 1, • ■ ■ , 2k. Up to an overall SO(2k ) rotation, this is the same as the exchange 
(2.37). Thus we have proved that the lambda parity induces the exchange of w m and v m . 


3 Tree amplitudes 


3.1 4-point 

The momentum conservation in the particle basis reads, 

P = —AiAi + A 2 A 2 — A 3 A 3 + A 4 A 4 = 0, (3-1) 

where we suppressed the spinor indices. In terms of the Lorentz scalars, (ij), we obtain 

(12) (23) (31) 


(34) (14) (24) 

(12 ) 2 + (23) 2 - (31) 2 = 0 


= 0 = ±1, 


(3.2) 

(3.3) 


The sign factor a in (3.2) specifies a branch of OG 4 . Without loss of generality, we will 
work in the <r = +1 branch for the rest of this subsection. 

In the light-cone basis, the gauge-fixed C-matrix and the metric are (u = m 2 ) 


C = 


10 0 u 
0 1-u0, 


9 = 


/o 0 1 o\ 
0 0 0 1 
10 0 0 
\o 1 0 0/ 


(3.4) 


In the particle basis, 


C = 


1 1 — u u 
u —u 1 1 


g = diag(-,+,-,+). 


(3.5) 


To avoid confusion, we put hats on the objects in the light-cone basis. 

In the evaluation of the Grassmannian integral, the kinematic delta-function gives 


S(C ■ A) = Jf 5 3 (P)8(u-m). 

The value of m is determined by (2.24): 

2{ Wl w 2 ) (23) - (31) (12) 


u* = 


(rciu 1 ) + (W2V 2 ) 


( 12 ) 


(23) + (31) 


(3.6) 


(3.7) 
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The equality of three expressions follow from (3.2) with a = +1 and (3.3). The Jacobian 
factor in (3.6) is 


Jf = ^ = _(( 23)+ (31)). (3.8) 

The fermionic delta function gives 

5(C-rj) = j[5 e (Q)= (^y)V(Q). (3-9) 

The denominator at u = u* is 


1 


AIi AI2 


1 _ ((23) + (31))- 

(-2«*)(l-«2) 4(12)(23) 


(3.10) 


Collecting all ingredients, we reproduce the standard form of the 4-point amplitude, 


5 3 {P)5 6 (Q) 
(12) (23) 


(3.11) 


3.2 6-point 


It is well known that the Grassmannian integral for the 6-point amplitude is fully localized 
by the delta functions and leaves no contour integral. In the particle basis, the gauge-fixed 
C matrix in the (+)-branch is 


C+ 


^ 1 1 -U12 U12 -U13 U13 

U\2 —U12 1 1 — U23 U23 

\Ui 3 U13 U 2 3 -U23 1 1 


g — diag(—, +, —, +, —, +). 


The kinematic delta-function can be transformed into 


(3.12) 


S(C • A) = J" 5 3 (P) 5 6 {u mn - u * mn ). 

The value of u* mn is determined by (2.24) and the Jacobian factor is Jq 
fermionic delta function gives 


1 


S(C ■ r,) = J£ 5*(Q)5 3 (( + ), C| = p 

The Jacobian factor from the fermionic delta function is 

■*='£) 3 ' 

with R defined in (2.24). The denominator at Uij = u*j is 


AJ+M+AT+ 


U mn—U rnn 


8(^13 U 12 U 23)( U 23 u l2 U l^)i U 12 U 13 M 23) 


(3.13) 
1/2. The 

(3.14) 

(3.15) 


(3.16) 
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Collecting all ingredients, we obtain the result of the Grassmannian integral in the u-gauge: 


_ 4g 3 (P)* 6 (QM 3 (C+) __ 

R 5 ( u 13 ~ U 12 U 23)i U 23 ~ u 12 u l‘s)(. u \2 ~ ^ 13 ^ 23 ) 

The full amplitude is the sum of the contributions from the two branches related to each 
other by A-parity. As explained in the previous section, for odd k, A-parity exchanges the 
two types of light-cone coordinates. In terms of the super-space variables, the exchange 
means 


9 m ^9 mI , I - = \e mnp u™ n 9 pI . 

The explicit form of the C-matrix in the conjugate branch is 



f -1 

1 

u 12 

u 12 

u 13 

n 13 

c~ = 

-u 12 

-u 12 

-1 

1 

u 23 

u 23 


y—u 13 

-u 13 

—u 23 

-u 23 

-1 

1 


Summing up the two terms, we obtain the full 6-point amplitude: 


M =b\P)b\Q)~ % 




R 3 \(n| 3 u 12 u 23)( U 23 U 12 U 13)( U 12 U 13 U 23 ) 


+ 


5 3 (C-) 


( ul 3 - u\ 2 u 23 )(u 23 - u\ 2 u\ 3 ){ul 2 - u\ 3 u 23 ) 


(3.18) 


(3.19) 


(3.20) 


The planar 6-point amplitude can be factorized by two 4-point amplitudes in three 
different channels. At first sight, it is not clear how the result (3.20) can exhibit the 
factorization properties. Remarkably, it is possible to show that the consecutive minors 
from the two branches combine to produce the desired physical poles. In the u-gauge, we 
have 

Mj + M, = ^^ P 123 5 M 2 = — ^ 4*234) ^ P 345 j (3.21) 

where p 2 j k ... = {Pi + Vj + Pk + ■ ■ ■ ) 2 . A proof of this relation and discussion on its gauge 
(in)dependence is presented in appendix A.l. 


3.3 8-point 
Contour integral 

As discussed in section 2.2, the general solution to C ■ A = 0 in a u-gauge is 

^mn(^) = ^mn P l^^mnpq'U^ ■ (3.22) 

The Grassmannian integral reduces to a contour integral in z through the relation 

S 8 (C • A) = j£ 5\P) J dz S 6 (u mn - u mn (z )), (3.23) 
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with Jg = 1/(2 R). Up to an overall sign, the full 8 -point amplitude is obtained when the 
contour separates the poles of M\ and M 3 from those of M 2 and M 4 [18]. 

The minors of C-matrix can be at most quartic in u mn (z). But, explicit computations 
show that all quartic terms can be absorbed into the square of the quadratic polynomial, 

'bmn(z)Upq(z') T U m p(z')’Uq n (z') T U m q(z')U n p(z') — Z €mnpq • (3.24) 

Similarly, all cubic terms can be rewritten as the same polynomial (3.24) times a linear 
combination of u mn (z). These two statements imply that all minors of C, including the 
consecutive ones, are quadratic in z: 

Mi(z) = aiz 2 + biz + Ci = cii(z - zf)(z — z^~). (3.25) 

The fermionic delta function produces 

3 3 

S 12 (c ■ n ) = S 6 (Q) JliAjz 2 + B IZ + Cj) , (3.26) 

where the fermion bilinears (Aj, Bj,Cj) are defined as follows: 

\e mnpq u* mn + \z€pr Xy ul y e rI ) {e{ + \ze qszw u™e sI ) 

— T f_^.2 j-.mnnplnql , mnpq -* -xyarl al , f mnpq,* fl/fl/l 

— 4 |_ * * l mnpq Ll * u u i ^prxy ^rnn Uj * u w q ' c u, mn u p u q J 

= A IZ 2 + B IZ + C! . (3.27) 


We wish to evaluate the contour integral 

dz Y\i = \{Aqz 2 + Bjz + Cj) 
c 27r?: nt=i ( a i z ‘ 2 + hz + Ci) 


m = f 


(3.28) 


Note that, as far as the evaluation of the integral is concerned, the distinction between 
(di,bi,Ci) and (Aj, Bj,Cj) is immaterial. To streamline notations, we replace (Aj, Bj,Cj) 
(I = 1,2,3) by (di,bi,Ci) (i = 5,6,7) in what follows. To illustrate the ideas behind the 
computation, we find it useful to introduce a family of integrals of similar form, 


f dz nj=n+2(«l^ 2 + bjz + Cj) 
J c 2 m n:^ 1 (diZ 2 + biZ + d) 


(3.29) 


We can recover the original integral (3.28) by choosing n = 3. 

These integrals share two crucial features. One is that they are homogeneous functions 
of the variables (aj,&j,Cj) with degree (—1) for * = l,...,n + 1 and (+1) for i = n + 
2,..., 2n + 1. The other is that they are invariant under the SL(2, C) transformation, 


ez + f 
gz + h' 



e SL( 2,C). 


(3.30) 
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It is instructive to consider the generators of SL( 2, C) one by one: 


Dilatation: z -> e t z , (3.31) 

Inversion: z—>—l/z, (3.32) 

Translation: z -> z + e . (3.33) 

The change in z can be reproduced exactly by the change in the coefficients: 

Dilatation: (a*, bi, d) —> (e*aj, bi, e~ f d ), (3.34) 

Inversion: (a*, bi, a) (cj, -bi, a*), (3.35) 

Translation: (a*, bi, Cj) —>• (a*, bi + 2aje, c* + 6*6 + a^e 2 ). (3.36) 


The integral (3.29) should be invariant under the SL(2, C) action on z, provided that the 
contour transforms accordingly. It follows that the result of the integral should be invariant 
under the change of coefficients listed above. 

To be specific, let us focus on the contribution of the contour C\ enclosing the two 
poles zf only. The residue theorem gives 


I n (C 1 ) = - ( ' Nn{z i +) 1 


where we defined 

2n+l 


<i\ \D n {z l ) D n (z l )/ z^ z ^ 

_ 1 _ N n (zf)D n (zi) - N n (zi )D n (zf) 

ai D n (zf)D n (zi) zf - zf 


n+1 


N n (z) = ]^[ Mj(z), D n {z) = JJ Mi(z), Mi(z) = diZ 2 + biZ + a . 
j=n +2 


i=2 


The product D n (zf)D n (z 1 ) is easy to evaluate. Using the relations 

+ . - h + _ ci 

*1+21=-> z i z i = — > 

Ol Ol 


we find 

D n (4)D n (zj-) = ar 2n n A 

where we defined short-hand notations 


n+l 


i=2 


11 


(3.37) 


(3.38) 


(3.39) 


(3.40) 


^ij — Cj bjd , At — d^j CjCti, 7 ij — dibj (ijb\ . (3.41) 

The new symbols (aij , f3ij , jij) obey simple SL( 2,C) transformation rules, 

Dilatation: (a^, fy, 7 ^) ->• (e -4 ^, Aj, e^-), (3.42) 

Inversion: (ay, Aj, 7y) ->• (7y> , (3.43) 

Translation: (a^-, Aj, 7*j) (a+? ~ 2 Aa + - 77ij) > (3.44) 
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so that A ij is fully invariant under SL( 2,C). 
To summarize what we have done so far, 


In (Cl) = 


nS An 


(«i) 


2 n-lNn( z l)Dn(z 1 ) N n (z 1 )D n (z 4 ) 


z; - z 1 


(3.45) 


The remaining ^-dependent part may look complicated as both N n (z) and D n (z ) are 
degree 2n polynomials in z. However, since we only need their values at the two solutions 
of M\ [z) = 0, we can take the polynomial quotients. If we denote the quotient and the 
remainder by 


N n (z) = Q n (z)Mi(z) + R n z + S n , D n (z) = Q n (z)M 1 (z) + R n z + S n , 


(3.46) 


the integral gives 


In (Cl) 


{ai) 2n ~ l (R n S n — S n R n ) 

nS An 


(3.47) 


The denominator nS* An has degree (2n; 2, • • • , 2; 0, • • • , 0). It remains to express the 
numerator F n = (ai ) 2n_l (R, n S n — S n R n ), which has degree (2n — 1; 1 - - - , 1; 1, - - - , 1), in 
an SL(2, C) invariant way. To do so, we introduce a few additional SL( 2, C)-invariants: 


Jijk — det( (i ;, bj , Cfc) — 6 o.[j 6 jC/ c j , 

R-ij — bibj 2 ( (ijCj -)- Cidj ), 

Lijkl — (Xij'Jkl ‘-ftijftkl T lijOLkl ■ (3.4b) 


For n = 0, the integral vanishes trivially as the contour can be pushed to infinity 
without encountering any poles. For n = 1, the numerator F n should be of degree (1; 1; 1) 
in three groups of variables and anti-symmetric with respect to the last two. It appears 
that J 123 is the only SL(2, C)-invariant with required properties. An explicit computation 
indeed shows that 


Fi — — Ji23 • (3.49) 

For n = 2, we look for a polynomial of degree (3; 1,1; 1,1) with total symmetry under 
permutations in the same group and anti-symmetry between the last two groups. The 
answer indeed respects all the desired properties: 

F'2 = \K 1{2 J m4 K 5)1 + ±K n Ji( 2 (4 ^ 3) 5) • (3-50) 

Finally, we turn to n = 3, our original problem. There are a number of ways to combine 
Kij, Lijki to construct SL(2, (C)-invariants with correct symmetry properties. Remarkably, 
the answer can be organized using only two such combinations: 

F 3 = —An Jl2 (5 Ji3 6 Jl4 7) + |A 1 (2 1 3J4)l( 5 A 1 6 1 7) • (3.51) 
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A remark is in order. The integrals I n (C) are defined in such a way that if we set, say, 
M 3 (z) = M^(z), I 2 (Ci) should reduce to I\{Ci). In terms of F n , we should have 

F 2 (12345)|“ 3 = 5 ” = A 13 Fi(124) . (3.52) 

The reduction does not look obvious from the expression (3.50). Similarly, it is not obvious 
how the reduction from / 3 to I 2 occurs: 

F 3 (1234567)|“ 4=7 ” = Ai 4 F 2 (12356) . (3.53) 

It is conceivable that the decompositions (3.50) and (3.51) are not unique, and some alter¬ 
native decomposition will make the reduction more obvious. 


8-point amplitude: the result 


In summary, the 8-point amplitude can be written as 


> 4.8 


^(P)J e (Q)(l + 7r)J 8 s jf 


F( 1) 


Al 2 Ai 3 Ai4 


+ 


F( 3) 


A 3 iA 32 A 3 4 


with the Jacobian factors 


•/> 


B 


1 

2R ’ 




(3.54) 


(3.55) 


With the A-parity operator ir, the (1 + 7 r) factor denotes the sum over two branches of the 
orthogonal Grassmannian. In (3.54), we removed the subscript from _F 3 and and added 
reference to the contour by F(i). Note that this form of the 8-point amplitude is valid in 
any u-type gauges. 


Physical and spurious poles 

The 8 -point amplitude (3.54) is the sum of two contour integrals, I(Ci) encircling the poles 
from M\ and J(C 3 ) encircling the poles from M 3 . Each term carries physical and spurious 
poles. The most convenient tool to analyze the pole structure is the on-shell diagram 
pioneered by [5] and elaborated for ABJM amplitudes in [23-25]. 


As = 




Figure 1. On-shell digrams for the 8 -point amplitude collect residues from Mi and M 3 . 
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Figure 2. Poles of I(C\) correspond to boundary components of the on-shell diagram. 

The on-shell diagram of ABJM amplitudes are built from two building blocks: a quartic 
vertex and an internal line. The quartic vertex is precisely the 4-point amplitude (3.11), 
and the internal line equates the kinetic variables of two legs from different vertices and 
integrates (/ d 2 l 3 A). 

The poles of the amplitude corresponds to boundaries of the on-shell diagrams. Each 
on-shell diagram has five vertices. Barring disconnected diagrams, each vertex yields ex¬ 
actly one boundary term. Figure 2 shows the five boundary terms from the on-shell diagram 
for /(Ci). Using the canonical coordinates for on-shell diagrams introduced in [23-25], we 
can easily see which consecutive minors vanish as we approach each of the five bound¬ 
ary components. To be specific, we adopt the coordinates of [24] associated with the OG 
tableaux. 




1 2 3 


Figure 3. OG tableaux for I{C\). 

The tableau for I(C \) is depicted in Figure 3. It can be translated to the C-matrix 
according to the rules explained in [24]. Let Cj be the z-th column of the C-matrix. We 
begin by setting the ‘source’ columns (Ci, C 2 , C 3 , C 5 ) to form an identity matrix. We assign 
a coordinate t v to each vertex. To fill in the ‘sink’ columns (C 4 , Cq, C 7 , Cs), we consider all 
paths from a source to a sink which may move upward and to the right but not downward 
or to the left. The path picks up ±sinh(C) if it passes through the vertex, or ±cosh(f„) 
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if it makes a turn at the vertex. The final matrix element is given by a polynomial of the 
form, schematically, 


a 


source,sink 


n f v - 

paths vertices 


(3.56) 


We refer the readers to [24] for details. All we need here is the remarkable fact that the 
consecutive minors are given by monomials of the sinh(t) factors.: 


Mi = 0, M 2 = S 1 S 2 , M 3 = S 2 S 3 S 4 , M 4 = S 4 S 5 , s v = sinh(t„). (3.57) 


In these coordinates, the boundary operation amounts to taking one of the coordinate 
variables to zero or infinity. The orientation of the untied diagram in the OG tableaux is 
shown in Figure 4. 



t (f —» 0) (f —» 00) 


Figure 4. Boundary operation in terms of canonical coordinates. 


To compare with the factorization channels in Figure 2, it is convenient to rescale the 
minors by an overall factor, M t = Mi/(s 2 S A ), 

Mi = 0 , M 2 = — , M 3 = s 3 , M 4 = — ■ (3.58) 

S4 S2 

The rescaled minor M 4 vanishes in the limit S 2 —;► 00 or S 5 —> 0. Through the prescriptions 
in Figure 3 and 4, the two limits give the two boundary diagrams on top of Figure 2, which 
in turn corresponds to the factorization channels for p 323 and respectively. 

In the contour integral obtained earlier, the simultaneous vanishing of Mi and M 4 , or 
equivalently the ‘collision’ of poles from Mi and M 4 , would result in the vanishing of A 14. 
It is then natural to expect that A 14 , a polynomial of kinematic variables, is proportional 
to p'f 23 ^ 567 - I n the u-gauge, we can explicitly verify the proportionality between A^ and 
physical poles. By symmetry, we expect that all of the eight physical poles are indeed 
associated with “collision” of roots of the minors: 


P123P567 OC -R 4A 14 , P234P678 « -R 4A 12 , p| 4 57 , 781 OC ^ A 32 , P^Pln « ^ A 34 • (3.59) 


The powers of R are fixed on dimensional ground. We leave the details of the verification, 
including the numerical coefficients, to appendix A.2. 

We can identify the poles for A 13 = A 31 in (3.54) as spurious poles. A standard 
argument in the Grassmannian integral uses the fact that 


A 8 = /(Ci) + I(C 3 ) = —/(C 2 ) - I(C A ). 
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Since A 13 = A 31 arises from I(C\) and /(C 3 ) but not from /(C 2 ) or /(C 4 ), it must be 
spurious. The physical poles (3.59), in contrast, appear in both contour prescriptions. A 
related observation is that the on-shell diagram for A 13 = 0 in Figure 2 can cancel against 
the same diagram from the boundary of /(C 3 ) if sign factors are properly assigned. 

We conclude this section with a few remarks on the generalization of the methods we 
used. The u-gauge has some advantages over more familiar gauges based on Euler angles. 
One of them is the decomposition of the fermionic delta-function, 


5 3k (C- v ) = 4 k 6\Q)H( I , J? k 


R 


c 7 = 


__ € rnn Pl ... Pk _ 2u * m j^ + ■ ■ ■ (% (k _ 2) + , (3.60) 


with u pq linear in the z coordinates in (2.33). Another advantage is that, as explained in 
appendix A. 2, the minors take a relatively simple form in the u-gauge. 

Finally, in anticipation of the generalization to 10-point or higher amplitudes, we note 
that the SL( 2,C) invariants are related to cross-ratios. For instance, 


Al4 _ (zi+ - Z 4 +)(zi- - Z4-)(zi+ - Z 4 -)(zi- ~ ^4+) 
KnK u (zi+ - zi-) 2 (z i+ - z 4 -) 2 


(3.61) 


Higher point amplitudes would inevitably give rise to more complicated “collision of poles” 
and it would be crucial to introduce higher dimensional analogs of A, J, K. L invariants to 
work without explicitly solving quadratic equations for the z coordinates. 


4 Soft theorem for ABJM amplitudes 

Soft theorems in gauge (gravity) theories explore the limit in which one or more gluon 
(graviton) approaches vanishing momenta. It is well-known that the soft limit of a non¬ 
vanishing tree amplitude is divergent and that the leading divergent term takes a universal 
form. More recently the sub-leading terms in the soft limit were calculated by using on-shell 
techniques [36] and spurred renewed interest in soft theorems and their applications. In this 
section, we derive a soft theorem for ABJM tree amplitudes, following a similar analysis 
for three-dimensional supergravity theories [32] (see also [23] for an early consideration of 
the double soft limit of ABJM theory). 

Since the ABJM amplitudes are well-defined only for even number of external particles, 
it is natural to define the double soft limit of the (2k + 2 )-point amplitude A 2 k +2 by scaling 
the momenta of the last two particles, 

( P2k+l,P2k+2 ) —t £ 2 (P2k+liP2k+2) , (4.1) 

and taking the e —> 0 limit. In spinor variables, the scaling rule is 

(A 2 fc+i, A 2 AH- 2 ) —> e(A 2 fc+i, A 2 /H- 2 ) ■ (4.2) 


- 18 - 





In view of the soft theorems in gauge theories in various dimensions, we anticipate that 
A 2k +2 in the soft limit reduces to the A 2 k up to a universal soft factor 5(e), 


-^■2fc+2 |e—>-0 — S(e)A 2 k . (4.3) 

We will find that the soft factor consists of a leading and a sub-leading term: 

5(e) = 4*S (0) + -5« . (4.4) 

e z e 

4.1 Recursion relation for soft limit 


Following the approach of ref. [32], we will use the BCFW recursion relation for ABJM 
amplitudes to analyze the double soft theorem. It is convenient to choose the two reference 
particle in the BCFW recursion to be (2k) and (2k+ 1), namely, neighboring soft and hard 
particles. The BCFW-shifted kinematic variables are given by 


A2 k = c\- 2 k + sA 2 fc+i, f} 2 k = crj 2 k + sr] 2 k+i , 

A2/C+1 = s\ 2 k + cA 2 fc + i , fj 2 k+i = srj2k + crfefc+i • (4.5) 


where c = coshf and s = sinht with c 2 — s 2 = 1 and z = c + s = e t . 

As explained in [32, 36] for soft graviton theorems, only one of the terms in the BCFW 
recursion formula contributes to the divergent soft factors. In our notation, the term is 
depicted in Figure 5. Let us briefly review why this is the case. The recursion formula 
schematically takes the form: 


k 

A 2 k+ 2 = ^ A 

1 =2 


(L)H(z+,z-) z(R) 

21 2 • /l 2(fc+2-0 1 


(4.6) 


where A and “4 2 (fe+ 2 _Z) are two sub-amplitudes evaluated at BCFW-deformed momenta, 
Pf is the momentum through the factorization channel, and 

H(z,y) = 5AA=I) (4.7) 

x z — y z 

is the BCFW kernel introduced in ref. [18]. 

When Z = 2 as in Figure 5, in the soft limit (p 2 k+i,P 2 k+ 2 ) -4 e 2 (p 2k+ i,P 2 k+ 2 ), Pf 
in Figure 5 becomes nearly equal to p\ such that 1 /p 2 ~ 1/p 2 diverges by the on-shell 
condition of p\ . For l > 2 with generic assignment of external momenta, 1 /p 2 remains 

hnite. The other factors A^/\ -4 2 (a-+ 2 _;), ^( z +i z -) a l so remain hnite. 

Focusing on the term in Figure 5, we need to evaluate the following amplitude 


H(z z ) 

Al 2 fc+ 2 |e —>0 = / d 3 r]f At(A 2 fc + i, A 2 fe +2 , Ai, kf)(z +)^——- A 2 k(Af, ■ ■ ■ , k 2 k)(z + ) 


Pl,2k+l,2k+2 


+ (z. 1 - <->• Z-) 

= T+) 

— • /i 2fc+2 ^ "^2fc+2 ' 


(4.8) 
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2k 



Figure 5. Factorization diagram 


In this diagram, the BCFW-deformed pf is given by 

Pf(z) =P2+P3-\ - \-P2k-l +P2k(z) = -Pl ~P2k+2 ~ P2k+l(z) . 

The on-shell values z + , Z- are the two solutions of the on-shell condition pf(z ) 2 = 
can solve this condition order by order in e by assuming that s = so + esi+e 2 S 2 + - 
c = Co + eci + e 2 C 2 + • ■ •, with c 2 — s 2 = 1. The solutions are given by [32] 

c± = 1 2 ~ a ± + ^( e4 ) > 

e 3 

s± = — ea± + — [(a± + a T )a± — (a± — a- F )/3 2 ] + 0 (e 5 ), 
e 2 

z± = 1 - ea± - —a± + 0(e 3 ) 


where a,j and /3j are defined by 

(1,2fc + 1) ± (1, 2k + 2} 


a± = 


, P± = 


(2k, 2k + 1) ± (2/a, 2fc + 2) 


<1,2 k) ’ ^ <1,2fc) 

As mentioned in [32], the (+) solution corresponds to a = +1, i.e. 

(l,f) = (2k + l,2k + 2), (f,2k + l) = (l,2k + 2), (/, 2k + 2) = — (2fc"+~l, 1), 

and the (—) solution corresponds to cr = — 1, i.e., 

<!,/) = -(2k + l,2k + 2), (f,2k + l) =-{1,2k+ 2), (f,2k + 2) = (2k + l,l). 


(4.9) 

0. We 
■ • and 


(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 


These relations will be useful when we calculate the soft limit. 

Note that, in (4.5), the first correction terms carry e 2 weight relative to the leading 
terms for bosonic variables, whereas the relative weight is e 1 for fermionic variables. In 
order to compute the leading and sub-leading terms of the soft limit, we need only the 
leading correction terms for bosonic variables, but we should keep track of leading and 
next to leading corrections for fermions. 


- 20 





We first focus on the A^ +2 contribution to (4.8). Upon expansion in powers of e, the 
propagator in the factorization channel gives 

1 


1 


1 


(pi + P- 2 k+i + V 2 k+ 2) 2 e 2 a+a-( 1 , 2k) 2 ’ 

and the BCFW kernel gives 

4( z - ~ !) 


H(z+,z-) = 


a- 




a + — ol- 

The bosonic factors in the 4-point amplitude reduce to 

(f,2k + l) = e(l 1 2k + 2), 

(/, 2k + 2) = -(2Ar+l, 1) = —e(l, 2k + 2 ), 

(/, 1) = (2k + 2, 2k + l) = e 2 (l, 2k + 2)/3+ . 

The fermionic delta function in the 4-point amplitude can be rewritten as 

r3 / (2k + l,2k + 2) (2k + l,l) ^ 

<U(<2) = d 3 77/ + ' mk+2 ~ hi 


(4.16) 


(4.17) 


(4.18) 


<2fe + l,/> 


<2fc + l,/) 


x 5 3 (-(f, 2k + l)fj 2 k+i + (/, + 2}i]2k+2 - (/, l)??i) 

= S 3 (r]f + x) x (—e 3 (l, 2k + 2) 3 ) x d 3 (ij 2 k +1 + ?? 2 fc +2 + e(-a + r] 2k + /3+hi)) • (4.19) 

The first factor absorbs the fermionic integral in (4.8). The remaining factors combine with 
the bosonic factors to produce 0 (1 ) and 0 (e 3 ) terms in the soft limit. 

Collecting all ingredients, we find the z+ contribution to the soft limit of A- 2 k+ 2 - 


A, 


(+) = £ 3 (h2fc+l +V2k+2) 

2k+2 2e 2 (l,2 k)a+P + 


A 2 k(Vl,V2, • > V2k-l,V2k) , 


where 


m 

. V2k+2 


1 e/3+ 
e/3+ 1 , 


m 

1 V2k+2 


m k 


1 —eaq_ 

-ea+ 1 , 


V2k 

,V2k+ 1, 


(4.20) 


(4.21) 


Expanding explicitly in powers of e, we obtain 




h 3 (4+i) 


(+) __ 

2 fc+ 2 2e 2 (l, 2 k)a+f3+ 

1 

+ 2e(l, 2k)a + /3 + 
where we introduced 


A 2k 

[\^ijk0 { + + b Z (0k+ 1 ) (/?+7?2fc+2,i - a + R 2 k+i,2k)\ A 2 k ■ 

(4.22) 


d 


£+ = -a + r] 2k + /3 + rji , Rij = r)i— . 

dVj 


(4.23) 
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2 k —1 


2 k- 1 




Figure 6. Soft limit from the on-shell diagram perspective. 


Note that R t j are Absymmetry generators of the ABJM theory. Note also that we could 
have obtained exactly the same result if we had chosen external particles (2k + 2) and 
(1) as the reference legs for the BCFW recursion. In this sense, the symmetry between 
(2 k, 2k + 1) and (2k + 2,1) has been restored. This is natural from the on-shell diagram 
perspective as illustrated in Figure 6. 

The computation of ^2k+ r z ™ (4-8) proceeds in the same way. The result is 


a: 


(-) 

2k+2 


5 3 (0 k+ i) 


~A2k 


2e 2 (l,2 k)a-pJ 

+ „ v- ~ \h eiJK ®k+^k+l£- + ^(@k+ 1) (/3—i?2fe+2,l + Ot-R,2k+l,2k)\ -^2fc • 


2e(l, 2fc)a_/3_ 


(4.24) 


with = — a_?? 2 fc + (3-V i- Note that A^A_ 2 and A 2 fc i 2 are related to each other by 
lambda-parity, h. 2 k +2 —> —h- 2 k+ 2 , in addition to the exchange (a+,/3+) -H- (a_,/3_). 

Finally, we add the two contributions to obtain the leading and the sub-leading soft 
factor of the double soft limit 


-WUo = (^s (0) + ) A 2 k + 0(1) 


where the leading and sub-leading soft factors are 


^°) = 


1 


S 3 (0 k+ 1) g(gfc+i) 

2(1,2/c) [ ct+/3+ a_/3_ 

1 


(4.25) 


(4.26) 


[^e/JA'#fc + l#fc + l£+ + <5 3 (#fc+i) (/3+^2fc+2,i - a+-R2fc+i,2fe)] 


2(l,2A:>a + /3 + 

+ „ V-7T“ [gC/JA^fc+l^fc+l^- + <^ 3 (^fc+l) (P-R-2k+2,l + a-^2fc+l,2fc)] ■ (4.27) 


2(1,2 k)a-/3- 

4.2 Soft limit of the 6-point amplitude 

For simplicity, we will use A-parity operator n which is given by 

7T : A 6 -4 -A 6 . 


(4.28) 
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If we use this operator, we can consider only one part of 6-point amplitude. 

The soft limit with particle 5 and 6 soft is realized in the light-cone basis as 

ws^-ews, v 3 —> ev 3 . (4.29) 

As we observed earlier, the bosonic kinematic invariants receive leading corrections at the 


e 2 order. So we can freely use the 4-point kinematic relations. For example, 

R = + {w 2 v 2 ) + ( w 3 v 3 )) = (12) + (34) + 0(e 2 ) = 2(12) + 0(e 2 ), (4.30) 

— 1 = u[ 2 u^ + u"[ 3 ul 3 + u 23 ul 3 = u\ 2 ul 2 + d(e 2 ). (4-31) 

In the soft limit, up to 0(e 2 ) terms, the minors become 

2u* 

Mi = 2 (u* 13 — u* l2 u 23 ) = 2u^ 2 (— u\ 2 u* 13 — u 23 ) = —M^WiV 2 ) + {w 3 w 2 )) 

= ^(-2f)({45) + (46)) = e (_!Mh^ /?+ . (4.32) 

9 y* 

M+ = -2(u 23 - « 12 «* 3 ) = 2u* 12 (ul 2 u* 23 + u* 13 ) = -^(-(u 1 ^) + (wiw 3 )) 

= ^(2e)«15> + (16)) = e (I<1|^) a+ , (4,33) 

M* = -21'.,;, - «J 3 i4j) = -2u; 2 . (4.34) 

Recall that the 6-point amplitude (3.20) contains two fermionic parts 

5 3 (C + ), 6\Q 6 ). (4.35) 

Neglecting 0(e 2 ) terms, we observe that 

<5 3 (C+) = K 2 ) 3 5 3 (o 3 + = (u* 12 ) 3 5 3 (6 3 + e(-a+ t? 4 + (3+m)) ■ (4-36) 

V U V2 U V2 J 

The second identity follows from 


-V~°i - = -u* 23 u\ 2 0i + u* 13 u 12 9 2 = ((wsv 1 )^! + (w 3 v 2 )0 2 ) 

ll 12 U 12 H 

= — 2(l2) ^ W3 I(I 2 ) ~ l 1 ))( ? ?2 + r?i) + («7 3 |(|4) - |3))(?74 + %)] 

= — 2(12) Id 2 ) - l 1 ^ 2 + l 4 ^ 3 _ 

= e(-«+^4 + P+Vi ), (4-37) 

where we used (super)-momentum conservation 

-11)7?! + 12)772 - |3) t ? 3 + 14)7?4 = O ( e ), (14) = (23) + 0 ( e 2 ), 

and Schouten identities. 
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We now move on to the 4 6 (Qe) factor. To check our result, it is better to start with 
our conjecture. From our recursion relation result, A 4 part gives the super-momentum 
conservation like 


Qi = + e/3+%) + 12)772 - 1 3)773 + 1 4)(774 - ea+775) 


= -|l> 
= Qa — 


(>n + -M,) + |2)lfi . , 3)ib + 1 4) („. e M±p) , B ) 

«15)|4)i)5 + <46)|1)j? 6 + <16>|4)r? s + <45)|l)t,s) 


= Q4 - (((14)|5) + (45)11)) 775 - «14)|6) - (16)|4))? ?6 + (16)|4)t/ 5 + (45)1 1 ) 770 ) 

= Qe - e 11) (775 + T7 6 ) - e|^j-|4)(775 + t7 6 ) 


— Q 6 • 


(4.38) 


The last equality holds on the support of (4.36). So we can conclude that the six-point 
supermomentum conservation becomes the four-point supermomentum conservation with 
next-leading soft correction. Finally, our 6-point amplitude becomes 


-4,6|e-s>0 


n f 32 S 3 (P)S 6 (Q 6 )S(( + )) 

(1 + 7F) l 1&M+M+M+ ) 

( 5^(63 + e(-a +m + /3 +m ))5 3 (P)5 e (Q 4 ) \ 
{ 7Fj ^ 2e 2 (14)a + /3 + (12)(14) ) 

(^ S(0) + “ 5(1) ) A4 ’ 


(4.39) 


if we expand the second line of above equations in terms of e up to leading and sub-leading 
orders. 


4.3 Soft limit of the 8-point amplitude 

In this last subsection, we examine the soft limit of the 8-point amplitude we computed in 
section 3.3. In view of the computational complexity, we content ourselves with checking 
the leading order soft factor . 

8-point amplitude with u-cyclic gauge 

To take the double soft limit of the 8-point amplitude, we revisit the computation of 
section 3.3 with two slight changes. The first is that, to be specific, we work in the u-cyclic 
gauge. The C-matrix is given by 

< i 1 -U12 u\2 -*13 7113 
c _ 7712 -7712 1 1 -7723 7123 

7713 -7713 7123 -7723 1 1 

\77l4 —U14 7124 —7724 7734 “7734 


~U \4 71 i 4 \ 
-7724 7124 
-7134 U 34 

1 




»(*) 


(4.40) 


1 
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In this gauge, the fermionic delta function reduces to 


3 1 

5 l \C t?) = JI ^e mn n(6 I m + u mi {z)O lI W n + u nj {z)O jI W P + u pk {^ kI W q + u q i{z)6 11 ) 


i=i 


4S 6 (Q) H(A iZ 2 + B 1Z + Cl) , 4 


1=1 


(4.41) 


The fermionic bilinear coefficients are 


Ar — Ti mn 

— A tmnpq 


0 pI 0 qJ , B! = le^e^yu^u^ 1 ^ , Cj = \e mnp B4J T p 9 r q . 

(4.42) 


The second, more important change compared to section 3.3 is that, in order to expose 
the soft limit more clearly (more on this below), we use the contours C 2 and C 4 instead of 
Ci and C 3 . Of course the two choices are equal up to an overall sign. In the notations of 
section 3.3, the result is 


A 8 = 5 6 (P)5«(Q)(l + 'K)jgj£ 


F( 2) 


+ 


F(4) 


A21A23A24 A41A42A43 

The A-parity operator 7 r acts on Als as 


<443) 


7r : — y -A s . (4.44) 

The numerators F{2) and F( 4) are given by 

F( 2) = —A"22«/21^ 5 «/23 6 ^24 7 ^ + “ L 2 ^ 2 3 «^4)2(5A 2 6 2 7 ) > (4.45) 

F{ 4) = -A 44 J 4 1 (5 J43 6 J42 7) + ^ 4 ( 1 4 3 J 2 ) 4 ( 5 L 4 6 4 7) , (4.46) 

The physical poles are related to the A ij factors through 

2 10 2 10 

A2I = — ^4^234 4*678 > A23 = ~ - J ^ P 345 4*781 ■ 

2 10 2 10 

A 41 = - ^4 4*?23 4*567 , A43 = - p 2 8 12 • (4.47) 

The derivation of these relations is essentially the same as the one given in appendix A.2. 
The factor A 24 corresponds to spurious poles. 

When we consider the limit in which particles 7 and 8 become soft, divergent terms 
come from A 21 and A 23 . If we use the contours Ci and C 3 as in section 3.3, the two 
contributions are divided into two different on-shell diagrams. But, if we use the contours 
C 2 and C 4 , both contributions come from the residues of M 2 (z) and we can ignore the 
residues of M 4 (z). 
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Soft limit of 8-point amplitude 

In the lightcone coordinates, the the double soft limit of the 7 and 8 is realized by 

IV 4 -> eu) 4 , r 4 -> er 4 . (4.48) 

In the e —» 0 limit, -«* n4 and h ™ 4 are of order e. As we discussed earlier, kinematic invariants 
receive e 2 corrections, so we can freely use the kinematic relations of the 6 -point amplitude. 
For example, the identity (2.29) in the soft limit implies that 

1 + u\ 2 u]^ + u 23 u 23 + = 0 + 0(e 2 ). (4.49) 

To the leading order in e, the coefficients of minors Mj(z) are given by 
a\ = 4eu 42 u 34 , 61 = 4(u 42 u 42 + 1 ), ci = Aeu\ 2 u* u , 

a 2 = 2 e(-u 24 + h 23 u 34 + u 12 u \ A ), b 2 = 2 (u* 12 uf + u 23 u 12 + u * 13 + h 43 ), 
c 2 = 2e(—« 2 4 + U* 23 U * 34 + u* 12 u * u ), 
a 3 = 4 eu^uf , b 3 = ^{u 23 uf + 1 ), c 3 = 4 eui 4 U 2 3 , 

G 4 = 2 (—h 43 + u\ 2 uf ), 64 = 2 e(u\ 2 u\^ + u 23 u^ + ' 034 'u 23 + u^h 42 + u 24 + uj 4 ), 

C 4 = 2(—U 43 + 1142 ^ 23 ) • (4.50) 


If we focus on the leading order only, the supernromentum-conserving delta function 
of Ag trivially reduces to that of Ae- 


S 6 (Q 8 )\^ 0 = S 6 (Q 6 ). (4.51) 

The only non-vanshing contribution from the fermionic part in the numerator F( 2) is 

-F(2)|e-s>o = — (b 2 ) 3 C 4 a\ 2 a 23 a 3 a^ + (6 2 ) 3 «47i2723 C 5 C 6 C 7 . (4-52) 

One can easily check that fermionic bilinears 04 - 1 -/ = Aj and C 4+1 = Cj become 

A/U0 = -\e mnp uT n o pI e AI = -c£ x e AI , (4.53) 

CVUo = \e mnp u* m J T p 0{ = Ci X e{ . (4.54) 

The soft limit of Ay gives 

A 2 i = —ai 2 7 i 2 , A 23 = —a 2 37 23 , A 24 = (b 2 ) 2 a4C4 . (4.55) 


We observe that the following useful identities hold in the soft limit: 

04 = —M 3 , 
c 4 = -Af+ , 


a 12 = — 


8eM +(68) x 0|) 

R z R z 


712 = -8e 


(67) - (68) 

R 


x M 2 = —8 e/3- 


( 16 ) 


R 


M, 


2 ) 


a 23 = 


8e <ro + m> X + = 8m+ M +J 


R 


R 


723 = 8e <HL_M x M - = , 


(4.56) 

(4.57) 

(4.58) 

(4.59) 

(4.60) 

(4.61) 
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where here denote consecutive minors of C Y± contributing to Aq. 

Collecting all ingredients, we obtain the soft limit of the 8-point amplitude in the 
leading order 


As 1 e— >0 — (1 + 7 r) 


F( 2) 


R A A21A23A24 


_ (1 _|_ ?r ) I ~(^2) 3 C4«12Q23 Q 5 Q 6 Q 7 + (&2) 3 Q47l2723 C 5 C 6 C 7 

'R 4 J (^2) 2 04C4«12a237l2723 


= (1 + vr) 


2 5 6 2 

i ? 4 


1 


<24712723 


- 0 , 5 ( 16(27 + 


1 


C4Q!i2a23 


-C5C6C7 


<1+,r) (s<M) ^MrM^L 3 - + {<+) ” ( ' 


= -S W U 6 + 0(l/e). 


, 0 -H- 0 } 

(4.62) 


In the hnal step, we used the following non-trivial identity 


b 2 R 

(16) 


z, - rt / * -23 , * -12 1 * 1 -13\ 

(^ 12 ^* + ^23^* + ^13 + w * ) 

2 

((^l^ 3 ) - (w 1 ^) + (wiw 3 ) - (uV)) =8. 


(4.63) 
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Appendix 


A Closer look at factorization channels 


A.l 6-point amplitude 


From the Grassmannian integral point of view, the three factorization channels of the 6- 
point amplitude originate from the product of minors from the two branches, when the 
minors are computed from the C matrix satisfying C ■ A = 0. In the alternating signature 
metric g = diag(—I-1-b), the precise relation can be written as 

M$*Mw = K(-l) <+i+ %k- (A.l) 

The subscripts ( ijk ) denote a set of three columns from the C-matrix, and the superscript 
(±) distinguishes the two branches. The prefactor K may depend on the gauge choice but 
is independent of the choice of columns (ijk). 

Since the relation (A.l) is fully gauge-covariant, it is sufficient to prove it in any 
convenient gauge. Moving to another gauge alters the prefactor K, but the relation remains 
intact. We find it convenient to begin with a slightly unfamiliar “A-gauge”: 

( \a \a \a in \a \a \ 

~ A 1 A 2 a 3 A 4 a 5 A 6 ] /a o', 

±(35) (46) ±(51) (62) ±(13) (24 )J ' { J 

This gauge was used in deriving a dual Grassmannian formula for ABJM amplitudes [29]. 

We can verify (A.l) by explicit computations. Two well-known facts reduces the 
number of independent terms considerably. First, (ijk) and its complement (ijk) have 
minors that are equal up to an overall factor [23, 24], In the (k,k) signature, they satisfy 


M ijk M ijk M ijk M ijk 


(A.3) 


Second, the gauge (A.2) respects a permutation symmetry among odd labels and another 
permutation symmetry among even labels. So, without loss of generality, we may only 
check (ijk) = (246) and (ijk) = (456). 

The first one trivially gives 


M+ 6 M 2 - 46 = ((24) 2 ± (46) 2 ± (62) 2 ) 2 = (p 2 46 ) 2 . (A.4) 


The second one is slightly more involved. 

M 456 M 456 = ((^4)(45) ± (56)(62)) 2 - (46) 2 (13) 2 

= (24) 2 (45) 2 ± (56) 2 (62) 2 ± 2(24)(45)(56)(62) 

- (46) 2 ((24) 2 ± (46) 2 ± (62) 2 - (25) 2 - (45) 2 - (56) 2 ) 

= ((24) 2 ± (46) 2 ± (62) 2 )((45) 2 ± (56) 2 - (46) 2 ) 

± (25) 2 (46) 2 - ((24)(56) - (62)(45)) 2 
= — ^246^456 ■ (A-5) 
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In the second step, we used momentum conservation. Note that the alternating signature 
metric is reflected in the square of a partial sum of momenta as 


Pm = (-1 + (-1 y +k W 2 + (-i)* +i <K> 2 . (A.6) 


In the last step, we used Schouten identity to cancel the last two terms. To sum up, in the 
gauge (A. 2), we have proved the relation (A.l) with K\ = p 2 46 - 

Let us compute the gauge dependent factor K in (A.l) for the u-cyclic gauge used 
in the main text. Other u-type gauges can be treated similarly. There are two ways to 
compute K. We may either compute it directly in a u-gauge or find a gauge transformation 
between a u-gauge and the A-gauge (A. 2). In the latter approach, we look for GL( 3,C) 
matrices G± satisfying 


C =b _ /Oib/^fzb 

X — U U 1I 5 


where we recall the C-matrices in the u-cyclic gauge, 



( 

1 

1 - 

-Ul2 

U12 ~ 

-U13 

U13 

Ct = 

U\2 - 

■U 12 

1 

1 - 

U23 

U23 


\Ul3 ~ 

-Ul3 

U23 

-U23 

1 

1 


f 

-1 

1 

u 12 

ZZ 12 

zz 13 

ZZ 13 

c- = 


-ZZ 12 

-ZZ 12 

-1 

1 

zz 23 

zz 23 



-ZZ 13 

-ZZ 13 

-ZZ 23 

-zz 23 

-1 

1 


U=U* 


The matrices Cr turn out to take a simple form 




1 / =FA? + =fA£ + A? TAg + Ag 

2 \ (35) + (46) (51) + (62) (13) + (24) 


(A.7) 


(A.8) 


(A.9) 


To make contact with the u-cyclic gauge, we note that the first two rows of Cr are light- 
cone coordinates and the third row can be rewritten as 

((35) + (46), (51) + (62), (13) + (24)) = |(^ 3 - u?,u* 31 - uf,u* 12 - uf) . (A.10) 

Using this and the identity (2.29), we find det(G ± ), 

det(G + ) = {R/df [-(« 2 3 - uf)uf - (u* 31 - ZZ 31 )ZZ 31 - («* 2 - u^uf] 

= {R/ 4) 2 [1 + « 2 ) 2 + (ZZ 23 ) 2 + (u 31 ) 2 ] = (R/ 4) 2 (M 2 - 46 ) u , 
det(G ) = (-R/4) 2 [(^23 — u 23 )u 23 + (u 31 ~ ^) u 3 i + (' u i 2 ~ ^* 2 ) u i 2 ] 

= (R/ 4) 2 [1 + K 2 ) 2 + C«23) 2 + («3i) 2 ] = (i?/4) 2 (M+ 6 ) u . (A.11) 

The determinants allow us to determine K u up to a sign, 

z4e = K\ = det(G + )det(G~)K u = (R/4) 4 (M 246 M+ 6 ) U = (R/4) 4 K 2 u (p 2 246 ) • (A.12) 
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To fix the sign and make a cross-check, we compute M + M directly in the u-cyclic gauge. 
We begin by rewriting p 2A Q in the light-cone variables, 

4 2 v\m = 4 2 «24> 2 + (46) 2 + <62) 2 ) 

= {v\ + Wi,V 2 + IU2} 2 + (V2 + W 2 , V3 + W3) 2 + (v 3 + W3, V\ + W\) 2 
= {v\2 + ^12 + ( W\V2) + (V1W2)) 2 + (cyclic). (A. 13 ) 

We introduced the short-hand notations, v mn = (v m v n ), w mn = (w m w n ). Using the 
identity R(v m w n ) = —v mp w np , the above identity can be further rewritten as 

(4-R) 2 p 2 46 = (Rv 12 + Rwi 2 + U23^13 - V 13 W 23) 2 + (cyclic) 

= R 2 (v 12 + W 12) 2 + {.V 23 W 13 - V 13 W 23) 2 + (cyclic) 

= R 2 (vi 2 + vf 3 + V 23 ) + R~{w \ 2 + w i 3 + ^ 23 ) + 2R 4 

+ (V23W13 - V 13 W 2 3 ) 2 + ('V32VJ12 - V 12 W 32) 2 + (V 31 W 2 1 ~ V 21 W 31) 2 

= (R 2 + v \ 2 + vf 3 + v 23 )(R 2 + w \ 2 + wf 3 + w 23 ) 

+ R 4 - (^ 12^12 + V 13 W 13 + V23W23) 2 

= (R 2 + v \ 2 + v \ 3 + v 2 3 )(R 2 + w \ 2 + u> 2 3 + w 2 3 ) = i? 4 M 246 M+ 6 . (A.14) 

In the second line, all cross terms cancelled out when summed over cyclic permutation. 
In the third and fifth line, we used R 2 = V 12 W 12 + V 23 W 23 + V 3 \ w 3 \. which is the same as 
(2.29). To conclude, we have verified (A.l) that in the u-cyclic gauge, 

M ^k M ^k = (|) (-1 ) i+j+k P 2 jk . (A.15) 

A.2 8-point amplitude 


In this section, we work in the u-factorization gauge introduced in section 2.2: 


< 1 U, 2 4 13 '414 1 -U12 -U13 -v-iA 

—U\2 1 U23 U24 U12 1 —U23 —U 24 

-U 13 -U23 1 U 34 U 13 U23 1 -^34 

\— ^14 —U24 —U 34 1 U \4 U 2 4 U 34 1 J 


/ U j m.rt. — 


»(*) ■ 


(A.16) 


As we mentioned earlier, the minors are quadratic polynomials in z. For later convenience, 
we organize the coefficients of the polynomials with some extra shorthand notations, 

Ml = z 2 + (ill 2 + ^13 + '&14 + U 2 3 + U 2 4 + 1134) + 1 

= (1 + V)z 2 + 2 Uz + (l + W), 

M2 = z 2 — (u 23 + U34 + U24) + (u 2 2 + U13 + U14) — 1 

= (1 + r 2 - s 2 )z 2 - 2(p 2 - q 2 )z - (1 + m 2 - n 2 ), 

M3 = z 2 + (u 2 2 + U34) — (uf 3 + M 2 4 + ^14 + ^23) + 1 

= (1 + r 3 - s 3 )z 2 + 2(p 3 - q 3 )z + (1 + m 3 - n 3 ), 

M4 = z 2 — (u 2 2 + U13 + U23) + (u 2 4 + U 2 4 + U34) — 1 

= (1 + r 4 — S 4 )z 2 — 2 (p 4 — q 4 )z — (1 + m 4 — 714 ) . (A. 17) 
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We introduced 


and 


v = 2 £(5?)/ 


-I 


£■ 




-hi 


r 2 

= (ll 23 ) 2 + 

(77 34 ) 2 + 

(li 24 ) 2 , 

S2 

P2 

= u* 23 ul 4 - 

17^ 4 77 43 + 77 34 77 42 , 

Q2 

?77. 2 

= ( u 2 3 ) 2 + 

(t7 34 ) 2 + 

(^ 4 ) 2 , 

77 2 

1‘3 

= («i 2 ) 2 + 

(n. 34 ) 2 , 


«3 

P3 

>k —34 1 4= —12 

= «12< + «34 w * , 


Q3 

777-3 

= K 2) 2 + 

(^ 34 ) 2 , 


n 3 

r 4 

= (^ 2 ) 2 + 

(77 43 ) 2 + 

(u?r, 

S 4 

P4 

= 77 42 77 34 - 

77^ 3 77 24 + 77^ 3 77 44 , 

Q4 

777.4 

= K2) 2 + 

K 3 ) 2 + 

(n 23 ) 2 i 

774 


hi 


S2 = COl 2 ) 2 + (Uf ) 2 + (u 44 ) 2 


*7 7 


q 2 = u\ 2 u 


34 


i24 


/* i 2 
‘131 

-:24\2 


- + u 14 

* \2 
12 / 

43 ) 2 

23 

T “23 “ 

+ K 4 ) 2 

24\2 
' w v ) 

<?4 = U 44 U 23 - 77* ' 1 “ 13 


77 23 , 


g 3 = u{ 4 u- + u 23 

= Ks ) 2 

- frM\2 


14 

ui — u 


‘24 


U2 


+ 77 34 77 


* ^2 
14/ ! 

i 4 ) 2 4 

u 24 

2 3 ) 2 ^ 
34\2 
* ) ! 

,-712 


23 \ 2 

* ) i 

,-713 


— u 2i u. 


(A.18) 


(A.19) 


n 4 — (itu ) 2 + ( u k ) 2 + ( u 34 ) 2 • 

Note that rh + Sh = V, ph + qh = U, rrih + = W for each h = 2,3,4. The quantities 

(aij , Pij , 7 y) can be expressed in terms of (m^, rih,Ph, qh, fhi Sh) in a simple way: 


«i 2 = -4 [(m 2 + l)g 2 - n 2 p 2 \ , 712 = -4 [p 2 (r 2 + 1 ) - < 72 ^ 2 ] , 

P 12 = 2 [{m 2 + 1 )(r 2 + 1 ) - n 2 s 2 ] , 

«i 3 = 4 [(m 3 + l)g 3 - n 3 p 3 ] , 713 = 4 [p 3 s 3 - g 3 (r 3 + 1)] , 

/?13 = -2 [(m 3 + l)s 3 - n 3 (r 3 + 1)] , 

«14 = -4 [(m 4 + 1)?4 - n 4 p4] , 714 = -4 [p4(?h + 1) - 94S4] , 

/?14 = 2 [(m 4 + 1)(?’4 + 1) - n 4 s 4 ] . (A.20) 


As explained in the main text, the two physical poles for p\ 23 and p| 67 are captured by 
A 14 . We can verify this by an explicit computation. With a slight rewriting to facilitate 
comparison with the u-gauge, we have 

(I) P123=(!) ((12) 2 + (13) 2 + (23) 2 ) 

= (1 + iq) (1 + m.4) — 2 (t7^ 4 ?7 44 + U 24 !!^ + 77 34 77 34 ) 2 + 54774 
+ ((^i 2 + u\ 2 ) + 2(u^ 3 7( 43 — h 43 772 3 )) (u 24 77 44 — u\ A U 24 ) 

+ ((«J 3 + «* 3 ) + 2(t7 32 77 42 - ul 2 U* 32 )) {ufu\ A - U^U^) 

+ ((ii * 3 + ^ 23 ) + 2(t7 31 7721 — 77 21 77 31 )) (t7 34 77 2 4 — 77* 4 77 34 ) 

= (1 + r 4 )(l + 7774 ) - 2(t 4 ) 2 + 54774 
+ (xi 2 + 27 / 12)272 + (x 23 + 2y 23 )z 23 + (x 31 + 2y 3 i)z 3 i, (A.21) 

and, similarly, 

(I) p36?= (£) (< 56 > 2 +< 57 > 2 +< 67 > 2 ) 

= (1 + 7’4)(l + 7774) — 2(t 4 ) 2 + S4774 

+ (-372 + 2yi 2 )zi2 + (-X23 + 2y 23 )z 23 + (-a: 3 i + 2y 3 i)z 3 i . (A.22) 
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Combining the two, we have 

2 8 . 2 
^ 4 ^ 123^567 = [(! + r 4 )(l + m 4 ) - 2(f 4 ) 2 + s 4 n 4 + 2 (yi 2 2 ;i 2 + 2/23223 + 2 / 31231 )] 

- (aJl2^12 + 3-23223 + ^3l23l ) 2 ■ (A.23) 

On the other hand, we deduce from (A.20) that 

1 9 

-A 44 = [(m 4 + l)(r 4 + 1) - n 4 s 4 ] - 4 [(7774 + l)g 4 - n 4 p 4 ] [f 4 (r 4 + 1) - </ 4 S 4 ] . (A.24) 

A lengthy but straightforward computation shows that (A.23) and (A.24) are equal. 
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